function x=sit1(A,b)
% rezolvare sistem inferior triunghiular
% Intrari: 
%		A=matrice sistem
%	        b=vector termeni liberi
% Iesiri: 
%		  x=solutia sistemului
	[n,n] = size(A);
	x = zeros(n,1);
% se calculeaza termenul initial;
	x( 1 ) = b( 1 ) / A( 1, 1 );
% se aplica formula de recurenta
	for i = 2 : n
		x( i ) = ( b( i ) - A( i, 1 : i-1 ) * x( 1 : i-1 ) ) / A( i, i );
        endfor

endfunction
